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ABSTRACT 

Context. Previous works have demonstrated that the generation of secondary CMB anisotropies due to the molecular optical depth 
is likely too small to be observed. In this paper, we examine additional ways in which primordial chemistry and the dark ages might 
influence the CMB. 

Aims. We seek a detailed understanding of the formation of molecules in the postrecombination universe and their interactions with 
the CMB. We present a detailed and updated chemical network and an overview of the interactions of molecules with the CMB. 
Methods. We calculate the evolution of primordial chemistry in a homogeneous universe and determine the optical depth due to 
line absorption, photoionization and photodissociation, and estimate the resulting changes in the CMB temperature and its power 
spectrum. Corrections for stimulated and spontaneous emission are taken into account. 

Results. The most promising results are obtained for the negative hydrogen ion H~ and the HeH + molecule. The free-free process 
of H" yields a relative change in the CMB temperature of up to 2 x 10~", and leads to a frequency-dependent change in the power 
spectrum of the order 10~ 7 at 30 GHz. With a change of the order 10~ 10 in the power spectrum, our result for the bound-free process 
of H" is significantly below a previous suggestion. HeH + efficiently scatters CMB photons and smears out primordial fluctuations, 
leading to a change in the power spectrum of the order 10~ 8 . 

Conclusions. We demonstrate that primordial chemistry does not alter the CMB during the dark ages of the universe at the significance 
level of current CMB experiments. We determine and quantify the essential effects that may contribute to changes in the CMB and 
leave an imprint from the dark ages, thus constituting a potential probe of the early universe. 

Key words, molecular processes - atomic processes - cosmology: early universe - cosmic microwave background - observations - 
theory . 
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' 1. Introduction of the universe, in spite of the smal l mole cular abundances. In 

• j ,r,,.T,\ ■ ^ f act ' tne ear hy work of Mao li et alJ £[994) suggested that the 

The cosmic microwave background (CMB) is one of the most molecular ities could smear out CMB fluctuations on the 

powerful tools of high-precision cosmology, as it allows one scale of the hOTizon> and at ±e „ time create new secondary 

to determine the cosmological parameters, the power spectrum fluctuations due to the i nter action with the velocity fields which 

of initial fluctuations and various other quantities. It is thus ^ t in to . clouds in the dark This work had 

important to have a detailed theoretical understanding of all ±e intention to explain why no CMB fluctuations had bee n 

effects that have a potential influence on CMB measurements, observed at that time 
Following the WMAP 3 year and 5 year resu l ts dSpergel et. al.L 

|2007t | Hinshawetal| |200l | KomatsuetaU |2008j) that con- Since thm ^ there has been considerable progress both in the un . 

firmed our standard picture of cosmology, we are looking d e rst anding of chemical processes in the early universe and the 

forward to the precise measurement that will be performed with molecular abundances , as well as in the interaction of molecules 

Planck 1 in only a few years. The measurement of the electron wkh ±e CMB and the generation of spectral-spatial fluctua- 

scattenng optical depth allows one to constrain the effective tions wh0e recombination was ori ginally examined bv lPe^blel 

reiomzation redshift and yields indirect information about an | 1%8) and Zerdovich- Kurt and Sunyaev ( 1969), and im proved 

epoch that cannot yet be observed. Recalling that the cross in severa i follow - up wo rks (Matsuda, Sato and Take|l \m& 

sections of bound electrons can be larger by orders of magnitude |j ones and WvseL flggl [sisaki and Takahari \m£ based 

compared to the cross section of free electrons, the optical depth Qn analytic methodSj today ' s com p U t e rs allow a detailed 

due to molecules may provide information on the dark ages ^atment of me recom bination process based on a reaction 

— : network that takes into account hundreds of energy levels 

http://www.rssd.esa.int/index.php?project=planck 
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for H, He and He + , as in the w o rk of Seager, Sasselov and 
Scott (Seager. Sasselov and Scottl l2000h . A simplified code 
reproducing the results of this detailed calculation is given by 
Seager, Sasselov and Scot t (JSSS). In a recent series of papers, 
Swit zer and Hir ata (2007a b.c) considered the recombination of 
helium in great detail. Deviations of the CMB spectrum from 
a pure blackbody have also been considered in various works. 
Dubr ovichl (1 19751) considered the effect of hydrogen recombi- 
nation hnes. iRubino-Martin, Chluba and Sunvaevl (120061 120071) 
examined d istortions due to helium and hydrogen lines in more 
detail, and IChluba and Sunvaev (200 6)) examined distortions 
due to the two-photon process. ISunvaev and Zel'dovichl d 19721) 
describe the rich ph ysics involved in the re combination process 
in great detail, and Wong, Moss and Scottl d2007l) give a good 
overview of recent improvements and uncertainties regarding 
the recombination process. 



Here, the main focus is on the postrecombination universe 
and possible imprints in the CMB from this period. The for- 
mati on of H2 during the dark ages has already been discussed 
by ISaslaw and Zipovl d!967l) . A more detail e d trea tment of 
molecules has been perform ed by iPuv et all (1 19931), Sta ncil, 
Lepp and Dalgarno dSLD98l) and by Galli and Palla (IGP98I) . A 
useful collection of analytic formulae for estimating the abun- 
dance of various mole cules after recombinat i on wa s given by 
Anninos and Norman dAnninos and Normanl 1 19961). Rec ently, 
this problem was re-exami ned by Puy and Signore dPS07l) . and 
Hirata and Padmanabhan (HP06) examined H2 formation in 
more detail, taking into account the effects due to non-thermal 
photons. 

Regarding the interaction of molecules with the CMB, the 
effects of various m olecules du e to th eir optical depths have 
been considered by IDubrovichl d 1994ft . The enhancement of 
spectral-spatial fluctuations due to the luminescence effect, 
which is well-kno wn from sta r s in re flection nebulae, has also 
been discussed by IDubrovichl d 1997b . Observation al prospects 
for He rschel and ODIN have been discussed by iMaoli et aU 
d2005l). and the releva nce for Planck has been assessed by 
Dubrovic h et al.l (120071) . The effects discussed were based on 
smear-out of primary CMB fluctuations due to the molecular op- 
tical depth and the generation of secondary anisotropics due to 
scatte ring with proto-objects in the universe. iMaver and Duschll 
(2005) provided a recent overview of different processes that 
may contribute to the opacity in primordial gas, and derived 
Rosseland and Planck mean opacities. 



Black] (120061) recently considered the influence of the bound- 
free transition of the negative hydrogen ion on the CMB, and 
found an optical depth of more than 10~ 5 at 10 cm -1 , which 
would have interesting implications for Planck and other CMB 
experiments. This contributed to our original motivation to 
examine this and other effects in more detail. As we will show 
below, however, the optical depth is due to H~ bound -free 
transitions is much smaller than reported by Black (2006), and 
must in addition be corrected for stimulated and spontaneous 
emission. On the other hand, there are other effects from H" and 
HeH + that are close to observational relevance. 



From a numeric al point of view, the pioneering work of 
lAnninos et al.l d 1997b provided a flexible and easily extendible 
scheme that is still widely used in state-of-the art simulations of 



the early universe, and which is also adopted in the public ver- 
sion of the Enzo co de (lO'Shea et all l2004t iBrvan et al.L 119951: 
Norman et al., I2007I) 2 . This scheme is extended here to account 
for effects during recombination and the evolution of primor- 
dial chemistry in the homogeneous universe. We further provide 
an extended overview of the essential processes that may influ- 
ence the CMB, determine the contributions from the most rele- 
vant species, and discuss the possibilities for detection with the 
Planck satellite. In section [3] we present the general picture re- 
garding the formation of the first molecules in the universe and 
give some analytic estimates for the abundances. In section|2] we 
provide a general discussion of the potential effects of primor- 
dial chemistry on the CMB. In section[4] we present the chemi- 
cal network for our calculation, which includes some new rates 
for the HeH + molecules that are given in the appendix. Section 
[5] explains the numerical algorithm for the chemical network, 
and section [6] presents the model abundances as a function of 
redshift. In section [7] we explain our treatment of the different 
species and discuss the observational implications. Further dis- 
cussion and outlook is provided in section[8] 



2. Imprints from primordial molecules on the CMB 

The interaction with the CMB is crucial to determine the evolu- 
tion and abundance of primordial molecules, and conversely, this 
interaction may also leave various imprints on the CMB photons 
while they travel through the dark ages. The most obvious im- 
print on the CMB is probably a frequency-dependent change in 
the observed CMB spectrum I(v) due to absorption by a species 
M. An upper limit on this effect is given by 



I(v) = B{v)e 



-tm(v) 



(1) 



where B(v) denotes the unaffected CMB spectrum and the opti- 
cal depth tm(v) of species M at an observed frequency v is given 
by an integration over redshift as 



v) = J dlcr M [v (l +z)]n M 

I fM{Z)(Tu [Vo(l +Z)] 

Jo 



n H ,oc 
Ho 



d+zf 



(2) 



-,dz, 



where <tm(v) is the absorption cross section of the considered 
species as a function of frequency, iim the number density of the 
species, n^.o the comoving hydrogen number density, c the speed 
of light, Ho Hubble's constant, /m the fractional abundance of 
the species relative to hydrogen, and z / the redshift at which it 
starts to form efficiently. Such an optical depth can be provided 
by resonant line transitions from molecules with high dipole mo- 
ments like HeH + and HD + , free-free processes or photodestruc- 
tion of species like H , He - or HeH + , which have a relatively 
low photodissociation threshold. However, Eq. Q] gives only an 
upper limit because absorption can be balanced by inverse pro- 
cesses (spontaneous and stimulated emission). A better estimate 
can be obtained by introducing an excitation temperature T ex de- 
fined by 



— = — exp 

m gi 



E„ - Ei 

kT e x 



(3) 



where n u and «/ denote the population of the upper and the 
lower level, g u and gi are the corresponding statistical weights 



http://lca.ucsd.edu/portal/software/enzo 
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and k is Boltzmann's constant. The excitation temperature is es- 
sentially determined by the ratio of collisional and radiative de- 
excitations. For molecules with non-vanishing dipole moments, 
the de-excitation is dominated by radiative transitions, and the 
level populations are in equilibrium with the radiation tempera- 
ture. The frequency-dependent change in the radiation tempera- 
ture is given by 



Continuum 



AT r ~ - (7V - T ex ) t m (v). 



(4) 



This expression indicates a major obstacle for the detection of 
primordial molecules: those species with high dipole moments 
and thus high cross sections have an excitation temperature 
which is very close to the temperature of the CMB, while 
species with low dipole moments have a very low optical depth. 
Molecular resonant line transitions are thus unlikely to lead to 
a non-negligible net change in the radiation temperature. For 
photodestruction processes, this is different because the destruc- 
tion is regulated by the CMB temperature, while the inverse 
formation processes are governed by the gas temperature. Thus, 
photodestruction can lead to a net change in the number of 
CMB photons. The same is true for free-free processes, which 
emit a blackbody spectrum according to the gas temperature, 
but absorb the spectrum of the radiation field. 



In addition, as discussed by [Maoli et al.| ( 1 19941) . the opti- 
cal depth from primordial molecules may smear out primary 
fluctuations in the CMB if the optical depth acts in a way that 
effectively scatters the CMB photons. This is possible even in 
a situation where stimulated and spontaneous emission balance 
the absorption of molecules, such that there is no net change 
in the number of photons. It is thus somewhat complementary 
to the effect discussed above. However, it must be noted that it 
requires spontaneous emission to be important, as stimulated 
emission does not change the direction of the photons, and does 
not provide a mechanism for scattering. We thus emphasize the 
importance to correct for stimulated emiss i on. A s shown by 
Basu, Hernandez-Monteagudo an d Sunvaevl d2004l) for small 
angular scales, such an optical depth (corrected for stimulated 
emission) then leads to a change in the power spectrum given by 



AC, ~ -2rC/. 



(5) 



where the C/'s are the usual expansion co efficients for the ob - 
served power spectrum. As discussed by Ma oli et alj (|l994), 
such scattering processes can in addition generate secondary 
anisotropies which are proportional to the optical depth of the 
scattering processes. However, this effect is suppressed by more 
than three orders of magnitude, as it is also proporti onal to the 
ratio o f the peculiar velocity to the speed of light. Dubrovich 
discussed the luminescence effect of various molecules 
which could potentially amplify the generation of secondary 
anisotropies. Based on the new abundances found in this work, 
we will give a basic estimate of this effect in section 181 iBasul 
(12007b considered in addition the effects of emission from 
molecules like HD and LiH + , as HD in particular is an important 
coolant in cold primordial gas. Unfortunately, the effect seems to 
be negligible. 

3. Recombination and the formation of molecules in 
the early universe 

3.1. Hydrogen recombination 

Various approaches exist to calculate the time evolution of the 
ionized hydrogen fraction during recombination. The traditional 



R 
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2S 



2P 



two-photon-decay 



Lyman-Alpha 



Ground state 

Fig. 1. The idealized three-level hydrogen atom (ground state 1 s, 
excited states 2s and 2p, continuum), and the relevant transi- 
tions. 



approach is based on solving one single ordinary differential 
equation for the ionization degree x e , which can be done using 
the appro ximation of an effective three-level-atom. This was first 
done by Peebles (1968), and an improved version of this equa- 
tion was derived bv lJones and Wvsa ([1985). This treatment was 
based essentially on the following assumptions: 

- It is sufficient to take into account only hydrogen (helium is 
not considered). 

- Collisional processes are outweighted by radiative processes 
and may be ignored. 

- The relative populations of the excited fine structure states 
are thermal. 

- Recombination and photoionization rates to and from higher 
states are related by Saha's formula. 

- Any such recombination cascades down to settle in the first 
excited state. 

- The level populations obey n\ s <k «2j, with n\ s the number 
density of hydrogen atoms in the ground state and «2j the 
corresponding density for the 2s state. 

- Each net recombination gives rise to a Lyman-ff photon or 
two others of lower energy. 



Seag er. Sasselov and Scottl (]2000) presented a very detailed 
calculation independent of these a-priori assumptions, treating 
helium and hydrogen as multi-level atoms with several hundred 
levels and evolving one ordinary differential equation (ODE) 
for each level, using a self-consistent treatment of the radiation 
field and its interaction with matter, which effectively leads to 
a speed-up of recombination. To allow the integration of their 
new method in other cosmolo gical applications like CMBFAST 
dSeliak and Zaldarriagal 1 1996b . ISSSl produced a simplified ver- 
sion of this code which is capable of reproducing the results of 
the full multilevel treatment. This simplified code, RECFAST 3 , 
solves three ODEs for the ionized hydrogen fraction x p , the ion- 
ized helium fraction xe e and the gas temperature T. These ODEs 
are parametrized so as to reproduce the speed-up in the recom- 
bination process. In this paper, we make use of the ODE solving 
for the ionized hydrogen fraction, which is given by 



[x e x p « H ffH -faO- - x p )e- h ^' kT ] 



dz H(z)(l + z)[l + K H (A H +/?h)«h(1 - * P )] 
x [1 + £hAh«h(1 - x p )l 



(6) 



http://www.astro.ubc.ca/people/scott/recfast.html 



4 



Dominik R. G. Schleicher et al.: Effects of primordial chemistry on the cosmic microwave background 



In this equation, «e is the number density of hydrogen atoms and 
ions, h p Planck's constant and the parametrized case B recombi- 
nation coefficient for atomic hydrogen an is given by 



an — Fx 1(T 



at" 



1 +cf d 



3 -1 

cnr s 



(7) 



with a = 4.309, b = -0.61 66, c = 0.6703, d = .5300 and 
t = T/10 4 K, which i s a fit by lPequignot et al.l (Il991h to the co- 
efficient of Hummer] (Il994l) . This coefficient takes into account 
that direct recombination into the ground state does not lead to a 
net increase of neutral hydrogen atoms, since the photon emitted 
in the recombination process can ionize other hydrogen atoms in 
the neighbourhood. 

The fudge factor F = 1.14 serves to speed up 
recombination and is determined from comparison with 
the multilevel-code. The photoionization coefficient /3h is 
calculated from the recombination coefficient as fin = 
au(27Tm e kT/hJ,) 3 ^ 2 exp(-h p vii,2s/kT). The wavelength Au^p cor- 
responds to the Lyman-Q' transition from the 2p state to the Is 
state of the hydrogen atom. The frequency for the two-photon 
transition between the states 2s and Is is close to Lyman-ar and 
is thus approximated by vh ; 2 S = cjXw^p, where c is the speed 
of light (i.e., the same averaged wavelength is used). Finally, 
Ah = 8.22458 s" 1 is the t wo-ph oton rate for the transition 2s- 
ls accor ding to Goldman ( 1989), H(z) is the Hubble factor and 
A"h = ■tf i2p /[8nH(z)] the cosmological redshifting of Lyman a 
photons. 



3.2. H 2 chemistry 

Due to the expansion of the universe, not all of the free elec- 
trons will recombine with protons. Instead, a freeze-out occurs, 
as recombination becomes less efficient at lower densitie s. The 
freeze -out abundance of free electrons was fitted by iPeeblesI 
(119931) as 



«H+ 
«H 



1.2 x 10 



-5 Q 



1/2 



hSl b 



(8) 



for cosmologies with total mass density parameter Qo, baryonic 
density parameter Q. b and h is the Hubble constant in units of 
100 km/s/Mpc. There is no contribution from helium to the free 
electron fraction, as helium recombines very early. The free elec- 
tron fraction leads to the for mation of molecules and ion s like 
H~ and Hj . As dis cussed by Annino s~and Normanl (1 19961) and 
lAnninos et al. 1( 119971) . their abundances can be calculated by as- 
suming chemical equilibrium, as their formation and destruction 
timescales are much shorter than the Hubble time. This yields 
the approximate expressions 



«H_ 

"H 
«H£ 

«H 



2xio- 9 r - 88 ^ 

«H 
«H 



3xl0 -14 r 1.8^ 



(9) 

(10) 



The matter temperature can be calculated by assuming that at 
Z > 200 it is the same as the radiation temperature, owing to ef- 
ficient coupling of the temperatures through the Compton scat- 
tering of CMB photons off the residual free electrons, while at 
Z < 200, where this coupling becom es ineffective, it evolves 
as fo r a simple adiabatic expansi on dSunyaev and Zerdovichl 
1 19721; lAnninos and Normanl 1 19961) . So 



with T = 2.726 K is the CMB temperature at z = 
dFixsen & MatheA 120021) . At redshifts z > 100, H is efficiently 
photodissociated by the CMB and H2 is mainly formed by the 
process HJ +H — > H2+H + . Assuming thatHJ is formed most ef- 
ficiently at redshift zo = 300 without being photo-dissociated by 
the CMB and that the hydrogen mass fract ion is given by fu = 
0.76, one obtains for the H2 abundance (lAnninos and Normanl 
1996) 



^ 2 x 10-°^ ( 1 +Z0 )^. 
"h nll b 



(12) 



3.3. Deuterium chemistry 

The HD abundance is mainly determined by the deuteration of 
hydrogen molecules (i. e. H2 + D + — > HD + H + ). It is thus crucial 
to have the correct abundance of D + , which is essentially deter- 
mined by charge exchange with hydrogen atoms and ions, i.e. 
the processes D + H + -> D + + H and D + + H -> D + H + . As 
will be shown below, D + is very close to chemical equilibrium, 
which yields the abundance 



«D+ 
«D 



1.2 x 10~ 5 exp(-43 K/T) 



1/2 



(13) 



when expression © is used. The abundance of neutral deuterium 
atoms «d can be determined by assuming that deuterium is al- 
most fully neutral, i. e. «d ~ ClbfoPc/mD, where fo is the total 
mass fraction of deuterium, ni£, the mass of one deuterium atom, 
p c = 3//g/87rG the critical density and G is Newton's constant. 
It is clear from expression ( fT3l that the abundance of D + drops 
exponentially at low temperatures. Direct deuteration of molec- 
ular hydrogen can thus only occur at redshifts where the expo- 
nential term is still of order one, before the exponential fall-off 
becomes significant. We thus evaluate the relative abundance at 
redshift 90. For both the formation and the destruction process 
HD + H + — > H2 + D + , we estimate the rates with the simple ex- 
pressions of lGP98l Again, we emphasize that more detailed nu- 
merical calculations should use the rates given in the appendix. 
The abundance is then given at z — 90 as 



«HD 

n H 



=90 



f Q 

1.1 x 10- 7 feexp(421 K/r ;=90 )^ 



3/2 



(14) 



As there is no efficient destruction mechanism for HD at lower 
redshifts, the fractional abundance remains almost constant for 
z < 90. 

3.4. HeH + chemistry 

As we will show below in more detail, chemical equilibrium is 
also an excellent approximation to determine the abundance of 
HeH + . For the rates presented in IGP981 the process of stimu- 
lated radiative association of H + and He dominates over the non- 
stimulated rate. With the new rates presented in appendix[B] we 
find that both rates roughly coincide for gas temperatures greater 
than 10 K. Thus, for an analytic estimate, we approximate the 
combined formation rate through stimulated and non-stimulated 
radiative association by taking twice the rate for non-stimulated 
radiative association. The dominant destruction process at low 
redshifts is charge-exchange via HeH + + H — > He 
yields for chemical equilibrium 



Hj, which 



r = r (i + z), z>2oo, T = j^r^( l+ d 2 > z< 200 ' (U) n ^ ~ L76x lo^VeOWn^f^^) ex p(-^o^j 
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7.03 x 1(T 16 «h 



Q}' 2 
_ 

h&B 



-0.24 



o / ^ / f \ The m ain difference compared to the deuterium network of 

Qb~ 1 300 K/ 6XP \ 4000 K/ IGP98I is the detailed treatment of D". We have added reaction 

44, estimating its rate from reaction 43, thus also considering 
the contribution of H . Als o, reac ti ons 27 and 28 are calculated 
fro m the revise d rates of ISavinl (|2002|) . As the fit provided 
by ISavinl d2002l) for the rate of reaction 27 becomes negative 
for T < 2.5 K, it is set to zero at these temperatures 4 . For 
temperatures larger tha n 200 K, the revised set of deuterium 
rates of Galli and Palla (GP02) are used. For lower temperatures, 
however, some of the new rates sho w unph ysical divergences. In 
these cases, we use the rates from IGP98I when the temperature 
drops below 200 K. 



4. The chemical network 

In our chemical network, we have included the formation 
paths of H2 and HD. The complete list of rates is given in 
table IB. 11 and some new rates that are relevant for HeH + are 
discussed in more detail in the appendix. H2 can be formed 
by two main channels, via the reactions H + H — » H2 + e 
and Hj + H — > H2 + H '*" . A very good com pilation for the H2 
chemistry was given by Yoshidaet al. (2006). Our compilation 
for the H2 formation rates is similar, but we do not include 
all of their three-body-processes, as they are not relevant in 
the low-density regime explored here. Also, we keep those 
modificati ons for low temperatures that were originally given by 
IAAZN971 The ionized fraction of hydrogen is not deter mined 
by solving rate equations, but from the RECFAST code of lSSSl 
The photodissociation rates for most of the hydrogen species are 
those of GP98. For the photodissociation of H^, we use the rate 
from their standard model, which assumes that the the levels of 
the molecule are populated according to LTE. Since the Hjlevel 
populations will be strongly coupled to the CMB at the redshifts 
at which H^photodissociation is significant, this assumption 
is more reasonable than using a rate that assumes that Hjis 
completely in the ground state. However, a better understanding 
of this molecule would certainly be desirable. For the molecule 
HD + , we estimate the photodissociation due to reactions 22 and 
23 of Table A.l as half of the corresponding reaction for Hi. 
This is in agreement with recent isoto pic helium experiments 
( Pede rsen et all 120051; iBuhr et all 120071) . which found a similar 
effect for dissociative excitation with helium isotopes. 

Photodissociation of H2 by the Solomon process (reaction 
20; Stecher and Williams, 1967) is calculated following the 
procedure described in IGJ07L with the assumption that the 
rotational and vibrational levels of H2 have their LTE level 
populations. 



In a recent work of Capitelli et al. dCCDL07l) . it was shown 
that the reaction rate for the process H2 + e — > H + H is 
larger by several orders of magnitude than the rate given by 
GP98J This is because they include the vibrational levels of 
molecular hydrogen in their calculation, and find that the 
excited vibrational levels cannot be neglected for this process. 
The change in the order of magnitude, however, does not lead 
to a significant change in the results, as this rate is multiplied 
by the rather small densities of molecular hydrogen and free 
electrons. 

The hydrogen and helium chemistry is almost completely 
decoupled, as can be seen from the small rates for the charge 
exchange reactions 50 and 5 1 . Their main interaction is via the 
HeH + molecule. Owing to its relatively high contribution to the 
optical depth for the CMB, we have included this molecule in 
our chemical network and present some new rate coefficients for 
it in the appendix. The molecule HD gives a contribution to the 
cooling and is mainly formed through by the deuteration reac- 
tions 30 and 32, but there are also contributions from reactions 
31, 43 and 44, involving HD + , D and H . For the formation 
of deuterated molecules, it is therefore important to determine 
the ionized fraction of deuterium, which is given through the 
charge exchange reactions 27 and 28. Our deuterium network i s 
inspired by the compilation of Nakamura and Umemura (2002). 



5. The numerical algorithm 

To determine the evolution of primordial molecules in the early 
univer se, we have employed the chemical netw o rk of the Enzo 
code dO' Shea et all 12004 iBrvan et all 119951; iNorman et all 
20071) and extended the numerical approach developed by 



Anni nos et all d 19971) for primordial chemistry to determine the 



chemical evolution of the homogeneous universe. The main is- 
sue was the calculation of the ionization fraction, which is de- 
termined by complex interactions between the CMB and the 
ground state as well as the excited states of atomic hydrogen, 
and which goes beyond typical applications of primordial chem- 
istry. We thus included the RECFAST code of lSSSl as a subrou- 
tine for this calculation. For the deuterium and helium species 
as well as the molecules, however, we use the fi rst order back- 
wards differencing (BDF) method developed by Annino s~et all 
d 19971) . The chemical timestep is set to 1% of the hydrodynam- 
ical timestep. The latter is given by Afh y di- = t7(Ajc/c s ), where 
Ax is the cell size, c s is the sound speed and 77 is a safety factor, 
here taken to be 0.5. This proved sufficient to resolve the rele- 
vant chemical timescales, and simulations performed with even 
smaller timesteps gave identical results. The rate equations for 
the species i are given in the form 



— = -DiHi + C„ 
at 



(16) 



where D,- and C, are the destruction and creation coefficient for 
species i, respectively, which in general depend on the number 
densities of the other species and on the radiation field. Equation 
[TSlis discretized and the right-hand-side is evaluated at the new 
timestep, yielding 



(17) 



At 

This can be solved for n" ew : 
C new Af + nf A 



(18) 

' 1 + D new Af 

The coefficients D" ew and C ( new are in general not known, but 
can be approximated using the species densities from the old 
timestep, and those spec ies from the n e w tim estep which have 
already been evaluated. lAnninos et all d 19971) has argued that 
H~ and Hj can even be evaluated assuming chemical equilib- 
rium, since their reactions rates are much faster than those of the 
other species. Assuming some species j in chemical equilibrium 
essentially means that fij = 0, yielding 

Cj 



n 



J,eq 



D; 



(19) 



4 Note that this does not introduce a significant error, as owing to its 
exponential dependence on temperature, rate 27 is tiny at T < 2.5 K. 
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Fig. 2. Results for the evolution of HJ and H~. 
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3. Results for the evolution of D + and HeH + 



In our chemical model, we have some additional species which 
have sufficiently fast reaction rates: HeH + , D and HD + . To 
check the validity of this assumption, we calculate the formation 
of molecules in two simulations: one using the non-equilibrium 
prescription ([T8T l for all species, and one using the equilibrium 
description ( fT9l for the species with fast reaction rates. The re- 
sults are presented in section [6] and confirm that the abundances 
of these species can be derived assuming chemical equilibrium. 



The BDF method is not a fully implicit numerical scheme, 
as several of the destruction and creation mechanisms must 
be approximated using the species densities from the previous 
timestep. However, we found that it is stable when the ionized 
fraction is provided from RECFAST. For consistency checks, 
we have varied the chemical timestep and explicitly ensured 
mass conservation for hydrogen, helium and deuterium, yield- 
ing consistent results. Note that mass conservation, charge 
conservation and po sitivity must be ensured as described by 
lAnninos et al.l d 1997b if the ionized fraction is not provided from 
an independent routine. 



We summarize 
way: 



- Loop over the chemical timestep t c i, em = 0.01f/ lw / ra until the 
species have been evolved through the total hydrodynamical 
timestep. 

- Update the temperature and the abundances of H, H + and e 
with RECFAST. 

- Calculate the abundances of the atomic and ionized helium 
species using the non-equilibrium prescription ([T8T i, and the 
abundance of the molecule HeH + using (fT9l . 

- Calculate the abundances of H and Hi using ( fT9] >, and the 
abundance of H2 using (fT~ST > . 

- Calculate D using CjJ]), D + , D and HD + using ([19), and 
finally HD using dH). 



1e-14 



1e-16 



.€■ 1e-18 



1e-20 



o 

CD 
CL 

CO 



1e-22 



1e-24 



1e-26 




Redshift 



Fig. 4. Results for the evolution of HD + and D 
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Fig. 5. Results for those species which freeze-out and are not in 
chemical equilibrium. 



6. Results from the molecular network 

The formation of molecules after recombination was calculated 
for a ACDM model with D. dm = 0.222, Q b = 0.044, C1 A = 0.734, 
H = 70.9km/s/Mpc, Y p = 0.242, [D/H] = 2.4 x 10~ 5 , where 
Qf/ m ,Qt,QA are the density parameters for dark and baryonic 



matter as well as dark energy, Hq is Hubble's constant, Y p is 
the mass fraction of helium with respect to the total baryonic 
mass and [D/H] is the mass fraction of deuterium relative 
to hydrogen. These are the parameters from the combined 
set of WMAP 3-year data and the d ata o f the Sloan Digital 
Sky Survey (SDSS) dSpergel et. all I2007I) . We expect only 
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Cosmology / Ref. 


[e/H] 


[H 2 /H] 


[HD/H] 


WMAP 3 + SDSS 


2.63 x 1(T 4 


4.23 x 10-' 


2.39 x 10~ lu 


(PS07) 


~ 10~ 4 


1.5 x 10~ 7 


~ 2x lO" 10 


(hpq6) 




~ 2 x 10~ 8 




Analytic approx. 


2.0 x 10~ 4 


2.9 x 10~ 7 


1.6 x 10- 10 


(GP981 


~ 2.4 x 10"* 


~ 1 x 10~ 7 


8x 10"" 



Table 1. Freeze-out of free electrons, H2 and HD for different 
cosmologies at z=100. Results of other authors and the analytic 
approximations are included for comparison. WMAP3 + SDSS 
refer to the cosmological paramete rs derived from the c o mbined 
set of WMAP 3 + SDSS data dSpergel et. all l2007h dPS07l : 
HP06;[SSS|) also refer to ACDM universes, while (GP98) refer 
to a flat CDM universe. 



Cosmology / Ref. 


[D + /H] 


[HeH + /H] 


WMAP 3 + SDSS 

(PSOT) 
Analytic approx. 
(GP98) 


8.02 x 10 lfa 
~ 1 x 10~ 19 
2.1 x 10" 20 
~ 5x 10~ 17 


9.98 x lO" 15 
4.6 x 10-' 4 
1.3 x 10~ 14 

~ 6.5 x 10~ 14 



Table 4. Freeze-out of D + and HeH + for different cosmologies at 
z=10. Results of other authors and our analytic approximations 
are included for comparison. WMAP3 + SDSS refers to the cos- 
mological para meters derived from th e combined set of W MAP 
3 + SDSS data dSpergel et. all 2007). (lPS07t iHPOfl ISSSh also 
refer to ACDM universes, while (GP98) refer to a flat CDM uni- 
verse. 



Cosmology / Ref. 


[e/H] 


[H 2 /H] 


[HD/H] 


WMAP 3 + SDSS 


1.92 x 10- 4 


7.97 x 10-' 


7.52 x 10- 1U 


(PS07) 


~ 6 x 10 5 


1.13 x 10-" 


3.67 x 10- 10 


(HP06) 




~ 6 x 10~ 7 




(SSS) 


1.92 x 10- 4 






Analytic approx. 


2.0 x 10 4 


2.9 x 10~ 7 


3.0 x 10 10 


(GP98) 


3.02 x 10~ 4 


1.1 x 10~ 6 


1.21 x IQ- 9 



Table 2. Freeze-out of free electrons, H2 and HD for differ- 
ent cosmologies at z=10. Results of other authors and the an- 
alytic approximations are included for comparison. WMAP3 + 
SDSS refers to the cosmological paramet ers derived from th e 
combined set of WMAP 3 + SDSS data dSpergel et. all |2007l) . 
dPS07HHP06t lSSSh also refer to ACDM universes, while (lGP98h 
refer to a flat CDM universe. 



Cosmology / Ref. 


[D + /H] 


[HeH + /H] 


WMAP 3 + SDSS 

(PS07) 
Analytic approx. 
(GP98) 


5.15 x 10~ 9 
~ 1.5 x 10-' 
3.7 x lO" 5 
~ 7.5 x 10-' 


4.27 x lO" 13 
~ 2.8 x 10- 15 
4.5 x 10~ 15 
~ 1.8 x 10~ 14 



Table 3. Freeze-out of D + and HeH + for different cosmologies at 
z=100. Results of other authors and our analytic approximations 
are included for comparison. WMAP3 + SDSS refers to the cos- 
mological para meters derived from th e combined set of W MAP 
3 + SDSS data dSpergel et. all l2007h . dPS07t IHP06I: ISSSl) also 
refer to ACDM universes, while (IGP98I) refer to a flat CDM uni- 
verse. 



marginal changes in the results for the recently published 
WMAP 5 -year data, as the cosmologica l parameters did not 
change significantly ( Komats u et all 120081) . We h ave performed 
the calcula tion for the HD cooling functions of IGP98L [FBPR 
and ILLA05L and find that the results are not sensitive to this 
choice. For our H2 cooling function, we use that of IGP98L but 
H2 cooling of the gas is never important, and our results would 
not cha nge significantly if we were to use the revised cooling 
rate of iGlover & Abell d2008l) . The detailed evolution of the 
various species is plotted in Figs. |2]|5] In Fig. [2] we give the 
results for HJ and H~, Fig. [3] shows the evolution of D + and 
HeH + , and Fig. [4] the evolution of HD + and D . For the other 
species, the results are given in Fig. [5] For several interesting 
species, we give the results at z = 100 and z = 10 in tables 
□H and compare them to the res ults of IGP98L |PS07l ISSSl and 
the analytical approximations of lAnninos and Norm an (1996) 



and section [3] As the analytic approximations do not take into 
account the H channel of H2 formation, they underestimate 
the abundance found in the numerical calculation by roughly a 
factor of 2 at redshift 10. 

The results clearly show the typical evolution of primor- 
dial chemistry as it is known from previous works. The two 
main channels for H2 formation, via H and Hj, are reflected 
in its cosmic formation history, yielding a first major increase 
at z ~ 300, where the relative abundance of HJ reaches a maxi- 
mum, and a second major increase at z ~ 100, at a maximum 
of the H fraction. At redshifts below 5, there is a new rise in 
the abundance of Hi. This is likely an unphysical feature from 
the fit to the rate, which is not valid below 1 K. However, the 
real evolution at these low redshifts will in any case depart 
from our calculation due to reionization, metal enrichment and 
structure formation, and this feature is not relevant with respect 
to the CMB. The evolution of the deuterium species essentially 
follows the evolution of the hydrogen species. Since deuterium 
and hydrogen are strongly coupled via charge-exchange re- 
actions, they recombine at almost the same time. However, 
due to the efficient charge-exchange reactions, there is no 
freeze-out of D + . Instead, its abundance drops exponentially 
at low redshifts. D and HD + peak at the same redshifts as 
H and HJ , and the evolution of HD resembles closely the 
evolution of HJ, as the dominant HD formation channel is given 
by H 2 + D + -> HD + H + . The evolution of HeH + consists of a 
first phase where its evolution is determined by the effectiveness 
of photodissociation, and a second phase where it is determined 
by charge-exchange with neutral hydrogen atoms. Due to the 
new formation rates presented in section [B] formation through 
stimulated and non-stimulated radiative association of He and 
H + is almost equally im portant and in total le ss effective than 
found in previous works dPS07HGP98 l :ISLD98h . 

Numerically, the free electron fraction found in this 
work roughly agree s with the analytic approximatio ns of 
lAnninos and Normanl (Qj96) and previous results of GP98, 
while |PS07| give a somewhat lower abundance. As the electrons 
act as catalysts for H2 formation, this is reflected also in the 
abundance of molecular hydrogen, and similarly in the abun- 
dance of HD , as it is primarily formed by direct deuteration 
of molecu lar hydrogen. At z = 100, the H2 abundance found 
by IHP06I is still more than one order of magnitude below the 
abundance we find here, due to the effects of non-thermal 
photons, while at redshift 10, this effect is much less important 
and their result is o nly 25% below the v alue found here. The 
analytic estimate of lAnninos et all d 19971) also underestimates 
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the H2 abundance at z - 10 by more than a factor of 2, as it does 
not take into account the H2 formation via H . At z = 100, the 
abundances for D + are of the same order of magnitude, although 
some differences exist, which may be due to the differences 
in the rates for the charge-exchange reactions, as well as 
differences in the abundance of ionized hydrogen. For the HeH + 
molecule, the abundance is still comparable to previous results 
at redshift 100, but lower by a factor of 5 at z — 10. 

We emphasize he re that the a bundance of H is similar to 
the one found by Black] d2006l) . The peaks in the abundance 
at z ~ 100 and z ~ 1400 are reproduced, and their height 
agrees with the results from our calculation when the physical 
number density given in his paper is converted to the fractional 
abundance given in Fig. [2] We estimate the uncertainty from 
reading off the values of his Fig. 1 to be a factor 2-3, and we 
point out that this is not sufficient to explain the difference in 
the optical depth that we find below. 

7. Effects of different species on the cosmic 
microwave background 

7.1. Molecular lines 

Due to the discreteness of the molecular lines, contributions to 
the optical depth arise only in narrow redshift intervals of the or- 
der Az/z ~ 10~ 5 , corresponding to the ratio between the thermal 
linewidth to the frequency of the transition. Peculiar motions 
of the order 300 x (1 + zY 1 ^ 2 kms -1 may further increase the 
effective linewidth by one or two orders of magnitude. This has 
no influence on the results presented here which depend only on 
the product of Az with the profile function that scales with the 
inverse of the linewidth, but may induce additional anisotropies, 
as we discuss in section [8] To evaluate Eq. [2] for a given 
frequency, we compute all the redshift intervals for which the 
photon frequency lies within an absorption line of the molecule, 
and add up the contributions from these frequencies. The 
relative importance of various molecules for the op tical depth 
calcul ation can be estimated with the formulae of iDubrovichl 
(11994 . The most promising candidate is the HeH + molecule, as 
it has a strong dipole moment and is formed from quite abundant 
species. Unfortunately, its destruction rate is very high as well. 
Another interesting molecule with a strong dipole moment is 
HD + . However, as a deuterated molecule, it has an even lower 
abundance. H2D + is no t considered because it has an even 
lower abundance (GP98), and has both a low abundance 
and no dipole moment. The other molecule in our chemical 
network with a non-zero dipole moment is HD. However, its 
dipole moment is eight orders of magnitude smaller than that 
of HeH + , and so even though its peak abundance is five orders 
of magnitude larger, its effects will still be negligible compared 
to those of HeH + . We therefore do not consider it further. In 
spite of its strong dipole moment, LiH is also not con side red, as 
it was already shown by Bougleux an d Gallil ( 1 19971) and [GP98 
that its abundance is lower by roughl y 10 orders o f mag nitude 
compared to the value assumed by iMaoli et all d 19941) . who 
discussed its potential relevance. 



For HeH + , we ca lculate the optical d epth by using the large 
dataset provided by lEngel et all d2005l) . which allows one to 
derive the line cross sections from the Einstein coefficients. 
The line cross section <tmj weighted by the level population 
for a transition from an initial state i with vibrational quantum 



number v", rotational quantum number J" and energy E", to a 
final state / with vibrational quantum number v' and rotational 
quantum number J' is given by 



o- M j(v) 



1.3271 x 10- 12 (2/' + l)c 2 



1 - exp 



QvrV 2 

h p v 
~kT r 



-exp 



AfMv-Vfd, 



-E" 



(20) 



Afi is the Einstein coefficient, v/t is the frequency at line centre 
of the transition i, 0(v - v/,) = AVf '^ ex P(~( v _ V/i) 2 /Av^) the 

profile function for the line width Avq - y/2kT/mMVfi/c, itim 
the mass of the molecule and Q vr is the partition function, given 
by Qvr = Zj gi ex P(~Ei/kT) with the degeneracies gi = 2J + 1. 
Again, we emphasize that the results are insensitive to the 
choice of the linewidth. The total cross section <tm is obtained 
as a sum over all cr M; . For temper atures between 5 00 and 
10000 K, we use the fit provided by lEngel etail d2005l) . while 
for lower temperatures, we do linearly interpolate between the 

values in their Table 5. The factor [l - exp(-h p v/kT,^ takes 
into account the correction for stimulated emission, which is 
especially relevant for the pure rotational transitions with low 
frequencies. As discussed in section [2] this must be taken into 
account as stimulated emission does not change the direction of 
the emitted photons. 



For HD + , we use the same formalism as for HeH + , but we 
determine t he partition function fr om the accu rate e nergy lev- 
els given bv lKarrandHilicol d2006l) . Foll owing IShuld 19911). we 
use th e transition moments | given bv lColbourn and Bunked 
d 19761) to determine the Einstein coefficients for the ro- 
vibrational lines as 



Afi 



32/rV},. 
3f\ p c i 



(21) 



where % p is the reduced Planck constant. 

The Einstein coefficients for the pure rotational transitions 
are calculated w ith the dipole moment Dq = 0.86 Debye of 
IDubrovichl d 1994 from 



32tt 3 i 



fi 



Di 



3 V 3 u 27+l' 



(22) 



7.2. The negative hydrogen Ion. 

There are two effects associated with the negative hydrogen ion 
than can affect the optical depth seen by CMB photons: the 
bound-fre e process of photodetachment that has also been dis- 
cussed by Black] (120061) . and free-free transitions that involve an 
intermediate state of excited H , i. e. 



H + e + y -> (FT)* H + e. 



(23) 



While the importance of the free-free process is well-known for 
stellar atmospheres, there has been little work on this proc ess in 
the low -te mperature regime. As the fit formulae given by iJohnl 
( 19880 and iGingerichl d 1961b diverge at low temperatur es, we 
have updated previous work oflDalgarno an d Land {1966) to cal- 
culate the free-free absorption coefficient for the low tempera- 
ture regime as described in appendix [Al and we use the fi t of 
IJohnl d 19881) to the calculation of lBell and Berringtonl d!987l) for 
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temperatures higher than 2000 K, where it is accurate with in 1 %. 
For the bound- free process, w e use the fit of lJohnl d!988h to the 
calculations of lWisharll dl979h . The treatment regarding absorp- 
tion, sponta neous and stimulated emission is based on the ex- 
pressions of iRuden. Glassgold and Sriul d 1990b . but with the up- 
dated cross sections mentioned above. Apart from the usual op- 
tical depth due to absorption, we introduce also effective optical 
depths due to stimulated and spontaneous emission. The follow- 
ing expressions have to be evaluated for the redshift z, and the 
frequency dependence is suppressed for simplicity. Of course, 
we emphasize that the observed frequencies at z — must be 
related correctly to the physical frequencies at higher redshifts. 
As usual, the contribution to absorption is given as 



dr, 



b f,abs 



+nft-(Tbfds, 



(24) 



where ds is the cosmological line element, <Tbf the cross section 
for the bound-free transition and «h- the number density of H . 
We further introduce the Planck spectrum Bj(v) of temperature 
T and frequency v. The contribution to the effective optical depth 
from stimulated emission is then given as 



dT bfMm = -(n H -) LTE o- bf e h " v/kT ds, 

where the LTE abundance (rtu-)LTE of H is given as 



A 



(nn-)LTE = n e n H —e' 



where A e 



(25) 



(26) 



^2nm e kT 

and vq = 0.754 eV the binding energy of H 



is the thermal de Broglie wavelength 
Note that the 



LTE abundance must be used here, as the processes of spon- 
taneous and stimulated emission depend on the actual density 
of electrons and hydrogen atoms, and in general the correc- 
tion for stimulated emission cannot be included as a factor of 
(1 - exp(-hv/kT)) unless t he H~ ion has its LTE abundance 
dRuden. Gl assgold and Shut Q990). The effective optical depth 
due to spontaneous emission is further given as 



dr 



bf, spon 



LTE 



2y 3 

c 2 B Tr (v) 



a- hf e 



-h 



? v/kT ds 



(27) 



For the free-free effect, stimulated emi ssion is alrea dy included 
in the rate coefficients a v (T) given by iJohnl dl988l) . which are 
normalized to the number density of neutral hydrogen atoms and 
the electron pressure. The effective contributions to the optical 
depth from absorption and emission are then given as 



dr 



'ff.abs 



— +nnn e kTa v (T)ds, 



dr 



ff.em 



— —nnn e kTa v (T)ds 



B Tg {v) 



(28) 
(29) 



The absorption by free-free transitions is thus proportional to 
a black-body spectrum for the radiation temperature, while the 
emission produces a spectrum determined by the gas tempera- 
ture. In the following, we will refer further to the optical depth 
from absorption, which we define as 



Tabs - I (dTbf,abs + dTffabs) , 



(30) 



HeHII photodissociation cross section 
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Fig. 6. HeH + photodissociation cross section obtained from de- 
tailed balance. 

which includes the corrections for emission. While the optical 
depth due to absorption, T flBS , is essentially responsible for pho- 
ton scattering and a change in the power spectrum according to 
Eq. (0, the effective optical depth r e /f leads to a net change in 
the radiation flux. The resulting change in the CMB temperature 
can be obtained by a linear expansion as 



AT = T ef f(v) 



B Tr (v) 



{dB Tr (y)ldT) Tr 



= f(v)Teff(v)T r , 



(32) 



where we have introduced a frequency dependent correction fac- 
tor f(v), which can be evaluated to first order as 



/(v) 



kT 

(l-expf-V/fcr])-^. 



(33) 



7.3. The negative helium ion. 



He~ is to some degree similar to the H . However, only the free- 
free process contributes in practice, as any bound states autoion- 
ize on a timescale of the o r der of hundreds of micro s econd s 
(Ho Menand MidtalL 1 1955b [Brage and Froese Fischerl 1199 lh . 
Thus, we take into account only the free-free process, which can 
be treated in the same way as for H~. We approximate the corre- 
sponding free-free coefficient by a po wer law prop ortional to v~ 2 
and normalize with the data given bv lJohnl d 1994b . Such a treat- 
ment should be sufficient up to frequencies of 1000 GHz, and 
thus for the frequency range interesting for the Planck satellite. 

7.4. Photodissociation of HeU + . 

We used detailed balance to determine the photodissociation 
cross section from inverse reaction. The latter is essentially de- 
termined by several narrow resonances, wh ich have been tab- 
ulated by Zygelman, Stancil and Dalgarno dZSD98h . The pho- 
todissociation cross section is thus given as the sum over reso- 
nances i as 



o-ph(v) = 2^ 



m e c (h p v — Eq) 



r,,,r ; /2 



(V) 2 



(h p v-E -E r j) 2 + (r/2) 2 



(34) 



and the effective optical depth 



T eff 



J {dT b j\abs 



+ dr, 



bf,stim 



+ dr, 



bf.spon 



where the parameters £V ; ,, F,- , and T, can be read off from Table 
2 of|ZSD98, and Eq is th e photodissociation threshold for HeH + . 
x From lDubrovichl 01997). we adopt the value Eq = 1.85 eV. The 
d T ff,abs + dTff tem J , (3 1 ) resulting photodissociation cross section is displayed in Fig. [6] 
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Fig. 7. The HeH + optical depth, both corrected and uncorrected 
for stimulated emission. The contributions from pure rotational 
and ro-vibrational transitions are given separately. 
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Fig. 8. The HD + optical depth, both corrected and uncorrected 
for stimulated emission. The contributions from pure rotational 
and ro-vibrational transitions are given separately. 



We consider only absorption, as it is sufficient to rule out the 
contribution from this molecule with respect to Planck. 
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Fig. 9. The absorption optical depth due to different processes: 
the free-free processes of H and He , the bound-free process of 
Ff~, and the photodissociation of HeH + . Clearly, the total optical 
depth is dominated by the processes involving H . 
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Fig. 10. The relative change in the CMB temperature due to the 
presence of H . We further plot the optical depths due to absorp- 
tion, spontaneous and stimulated emission for the bound-free 
process, as their overlap explains some features in the temper- 
ature change. 



7.5. Observational relevance and results 

As discussed in section [2] the optical depth from resonant 
scattering can effect the CMB by a change in the power spec- 
trum according to Eq. ©, and it may also produce secondary 
anisotropies. We neglect the latter effect for the moment, as it is 
suppressed by the ratio of the peculiar velocity over the speed of 
light. To quantify the importance of HeH + and HD + for a change 
in the power spectrum, we have calculated the optical depth and 
corrected for stimulated emission as described in the previous 
subsection. The results are given in Figs. [7] and [8] We find that 
the correction for stimulated emission is especially important 
for the lower freq uencies of the pure rotational transitions. 
As ex plained by Bas u. Hernandez-M onteagu do and Sunvaevl 
(2004), the sensitivity is not limited by cosmic variance when 
power spectra at different frequencies are compared, but the 
limit from instrumental noise corresponds to optical depths 
of 10~ 5 for the high-frequency bins. Thus, the signal is likely 



below the sensitivity of the Planck satellite by two orders of 
magnitude, but reasonable upper limits on the abundance of 
HeH + are feasible. From Figs. B.l.a and B.2.a in the appendix, 
we estimate the uncertainty in the formation rate to be a factor 
of 2, while the destruction rate might be larger by up to an 
orde r of m agnitude, if the old values of Roberge and Dalgarno 
dRDl |1982|) are adopted. This defines the main uncertainty in 
this result. Even with an improved instrumental sensitivity, very 
accurate foreground subtraction would be required and may 
create additional noise. 



Fig. [9] shows the optical depth due to absorption by several 
free-free and photodestruction processes. We find that it is 
dominated at low frequencies by the free-free contributions of 
H~, and at high frequencies by the bound-free process of H . 
The optical depth from the free-free processes of FT and He is 
essentially proportional to v~ 2 , at least for frequencies smaller 
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than 1000 GHz. This is essentially due to the characteristic 
frequency dependence in the absorption coefficient. Due to the 
approximation used for the Hefree-free absorption coefficient, 
the corresponding optical depth is somewhat overestimated for 
larger frequencies, but it is still only a subdominant contribution 
to the total optical depth. The effect from the bound-free 
transition d ominates, but is significantly lower than previously 
reported by iBlackl d2006l) . and the reason for this discrepancy 
is not obvious: as pointed out in section [6j our H~ abundance 
agrees with his within a factor of 2 or 3, whereas the optical 
depth at high frequencies is different by three orders of magni- 
tude. Although we have been unable to identify the reason for 
this disagreement, we have carefully checked our result and are 
confident that it is correct. 



As discussed in section [2] photodestruction and free-free 
processes can change the net number of CMB photons. The 
calculated change in the CMB temperature due to the free-free 
and bound-free effect is given in Fig. [10l and depends on the 
effective optical depth defined in equation (|3TI ). The change in 
the temperature due to the free-free transition is significantly 
lower than the absorption optical depth in Fig. [9] as the main 
contribution to the absorption optical depth comes from high 
redshifts z > 300, where the difference between radiation and 
gas temperature is very small and the optical depth due to 
absorption and emission balance each other. For the bound-free 
transition, the absorption is balanced to some extend by the 
spontaneous emission, and for frequencies larger than 1000 
GHz, spontaneous emission in fact dominates, as the CMB 
flux at these frequencies is low by comparison. Near 150 and 
1000 GHz, the contributions from spontaneous emission and 
absorption due to the bound-free process in fact become equal, 
making the net effect almost zero in a small frequency range. 



The net change in the CMB temperature due to the free-free 
and bound-free processes of H is small because absorption and 
emission processes are close to equilibrium. The optical depth 
due to absorption is thus considerably higher than the effective 
optical depth defined in Eq. (13T1 >. and can lead to a change in the 
power spectrum according to Eq. (0. Finally, we note that at red- 
shifts where reaction 2 of Table B. 1 dominates the destruction of 
H , there is an uncertainty of up to an order of magnitude in our 
predicted H abundance, owing t o the uncertainty in the rate o f 
reaction 2 discussed in detail in G lover. Savin & Jappsen (2006). 
However, this uncertainty does not appear to significantly affect 
the size of our predicted signal, as the dominant contribution 
comes from redshifts at which reaction 2 is unimportant. 

8. Discussion and outlook 

We have provided a detailed network for primordial chem- 
istry and solved for the evolution in the homogeneous uni- 
verse, and we examined the various ways in which primor- 
dial species can influence the CMB. The detailed calculation 
in the previous sections suggests that the H ion is only one 
order of magnitude below the detection threshold and strong 
upper limits on its abundance seem feasible, even though 
an accurate subtraction of frequency-dependent foregrounds 
will be required for this purpose. The relative deviations of 
the CMB from a p ure blackbody hav e been constrained by 
iMather et al1dl994 : lFixsen et alld 1996b to less than 1.5 x 10~ 5 . 
Distortions are also expected due to the two-photon process 



during recombination (Chluba and Sunyaev, 20 061) and the he- 
lium a nd hy drogen lines dRubino-Martm. C hluba an d Sunvaevi 
2006, 12007b . Finestructure transitions in heavy elements are 
also expect to produce some scattering af t er rec ombination 
dBasu. Hernandez-M onteagudo a nd Sunvaevi l2004t) . An accu- 
rate measurement of distortions in the CMB and a precise mea- 
surement of the CMB power spectrum can thus improve our un- 
derstanding of various processes during and after recombination 
if an accurate foreground subtraction is feasible. Given the tiny 
change in the CMB temperature found in this work, we further 
conclude that detecting the change in the power spectrum caused 
by scattering (see Eq. (|5]l) is the most promising way to obtain 
constraints on the chemistry of the dark ages by future CMB 
experiments. 

So far, we have only taken into account effects arising from 
the homogeneous universe. One might argue that the molecular 
abundances could be very different if there were a protogalaxy 
at the redshift of resonance. In fact, the redshift intervals for 
which resonance occurs are very narrow, corresponding to some 
100 pc. This might b e further increased by peculiar motions 
and lo cal turbulence. Nunez-Lopez, Lipovka a nd Avila-R eese 
d2006b demonstrated that HD in cold collapsed clouds can lead 
to a strong local fluctuation of the order 10~ 5 . Given the small 
volume fraction of such clouds, we neglect their impact in 
the present paper, although we will examine the effect of such 
inhomogeneous fluctuations in future work. 



iDubrovichl d 1997b suggested luminescence as an additional 
effect that may amplify the signal and lead to a strong frequency 
dependence. This effect is well-known for stars in reflection 
nebulae. It occurs at redshifts z = 300 - 100, where the rotational 
lines lie in the extreme Rayleigh-Jeans wing of the CMB and 
the first vibrational line near its maximum. Inelastic scattering 
in local velocity fields might thus provide much stronger 
frequency-dependent fluctuations. It is convenient to estimate 
the effect with the formula given by iDubrovichl d 1 997h : 



A7V 



A7V 



fm 



Vn 



10- 10 30km/s 0.1 



(35) 



where f m is the fractional abundance of the molecule, V p the 
peculiar velocity which can be estimated as V p = V p (0)/ yl + z, 
V P (P) = 600 km/s, and the qua ntity (AT r /T r ) n c an be conve- 
niently read off from Fig. 2 of IDubrovichl d 1997b . For HeH + , 
it yields a maximum effect of roughly AT r /T r ~ 10~ n , and 
AT r /T r ~ 10~ 12 for HD + , which are clearly below the sensitivity 
of the Planck satellite. 



As discussed by lLaunav. Le Dourneuf and Zeippenl d 199 lb . 
H2 is in general formed in excited states, but quickly decays 
into the ground state. Thus, H2 formation produces additional 
photons that may lie within the CMB radiation and thus produce 
distortions to the blackbody radiation. Prelimary estimates 
based on the transition between the first excited vibrational state 
and the ground state indicate that the effect yields a relative 
change in the CMB temperature of the order 10~ 15 , and is thus 
negligible. 

At the end of the dark ages, reionization will dramatically 
change the chemical evolution of the intergalactic gas, and 
produce large regions of ionized gas. In such regions, other 
formation channels for molecules could be relevant, like 
He + + H — » HeH + + y. The reaction rate of this channel is 
four orders of magnitude larger than the radiative association 
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of H + and He. On the other hand, the influence of destruction 
processes, such as dissociative recombination and photodis- 
sociation, will also be enhanced. The details of the evolution 
will also depend on the ratio of stellar to quasar sources and 
the details of the transition from Pop III to Pop II stars. A 
detailed analysis of this contribution is beyond the scope of 
this work, but clearly there is the possibility that this epoch 
could further increase the optical depths of HeH + and H . 
Other co ntributions may arise from heavy element s at low 
redshift. Basu, Hernandez-Monteagudo and Sunvaev (2004), 
iHernandez-Monteagudo. Verde and Jimenez! ( 20061) and iBasul 

(2007) suggested that one could use the change in the power 
spectrum of the CMB to constrain the chemical evolution 
of the low redshift universe. In fact, during the formation of 
metals and in early star forming regions, additional effects 
may occur tha t leave an interesting imprint i n the CMB. As 
described by Hernand ez-Monteagudo et ail (|2006), oxygen 
pumping may change the CMB temperature in metal enriched 
environments in a similar way as the Wouthu ysen-Field ef- 
fect t hat is well-k nown from 21 cm studies dWouthuvseni 
Il952t iFieldl 119581) . and the inhomogeneous distribution of 
metallicity in bubble-like struct ures may influence the CMB 
power sp ectrum as described by [H ernandez-Monteagudo et al. 

(2008) . IHernandez-Monteagudo. Rubino-Martm and Sunvaevl 
d2007l) further studied the effect of resonant scatter- 
ing during reionization and recombination. In addition, 
star-forming regions may perturb the primordial sig- 
nal through dust and molecular emission , espec ially CO 
(Rig hi. Hernandez-Monteagudo and Sunvaevl l2008bllal) . Such 
a potentially rich phenomenology will of course require a very 
careful analysis and a clear assessment of the different frequency 
dependences of various effects once the required sensitivity 
is reached. In the me an time, the increasin g sensitivity in 
instruments like Planck (Bersanel li et al.L fl99 6). the South Pole 
Telescope 5 (R uhl et al. , 2004) and the Atacama Cosmology 
Telescope 6 dFowler et al.L 120051) will allow at least to set upper 
limits that may constrain theories involving the dark ages, 
reionization and recombination. 
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Fig. A.l. The free-free absorption coefficient of H for 10 2 , 10 3 , 
10 4 and 10 s GH z, as a function of tem perature. Given are the fits 
of iJohnl (1 19881) and iGingerichl d 19611) for the high-te mperature 
regime, the calculation of IDalgarno and Land 0966) based on 
effective range theory as well as the new calculation of this work 
for the low-temperature regime up to 2000 K. 



(1966) used effective range theory (O'Mallev, Spruch and Rosenberg, 1961) 
to expand the cross sec tion, we use the more accur ate result of 
IDalgarno. Yan and Liul i 19991) as fitted bv lPinto and Galli |2008|), given by 



q (E) ■■ 



4x 10~ 15 cm 2 
(1 + £/(3.8eV))' 



(A.2) 



In the range between 1 and 10 7 GHz and between 0.1 and 2000 K, our results 
are well-fit by the expression 



Appendix A: Free-free transitions involving H 

While the importance of free-free transitions is well-known in stellar atmo- 
spheres, this process is usually not considered in colder environments. Thus, high 
precision calculations are only available for stel l ar temperatures b e tween 1000 
and 10000 K IJohnl 119881 ; iBell and Berringtonl , 119871 ; IGingerichl , 119611) . The 
free-free absorption coefficient kff is normalized per hydrogen atom and per unit 
electron press ure, and con tribut es to the optical depth by dr = kffnnn e kTds. 
As the fits by John 1 1988) and Gingerich 1 1961) diverge at low temperatures, 
thus giving an unphysical high contribution at low redshifts, we calculate the ab- 
sorption coefficient based on the formalism of Dalgarno and Lane 1 1966), which 
takes into account only the contribution of the leading term to the transition mo- 
ment. This approximation is valid at low energies and is thus a reasonable choice 
for the low-temperature regime. Assuming that the initial electron energies are 
described by a Maxwell distribution, the free-free coefficient kff is given as 

% (v) = 9.291 x 10-7-/ 2 (1-^)JTV (l + ¥)' /2 

+ q (E + h p v))dE cm 4 dyne -1 , (A.l) 

where v is the photon frequency, qo(E) the zero-order elastic scatteri ng cross 
section and a = h p v/kT a dimensionless parameter. While Dalgarno and Lane 



k ff (v) = lo( ni +^+^ 2 )(*i +*y+z> 3 v 2 ) cm 4 dyne -4 (A 3) 

where x = log] (v/GHz) and y = log[ (T/K), and a\ = -3.2421, a 2 = 
-0.502052, a 3 = 0.0117164, b x = 4.05293, b 2 = 0.169299, 63 = -0.00548517. 
At ste llar temperatures , our r esults differ by about 10% from the calculation 
of Bell and Berrington 11987), which is based on a more detailed treatment, 
whereas at lower temperatures, we expect an even higher accuracy of our result. 
We thus adopt the fit of John ( 1988) to the calculation of Bell and Berrington 
i 19871) for T > 2000 K and equation )A.3t for lower temperatures. 



Appendix B: Reaction rates 

Regarding the formation and destruction of HeH + , we have updated some of the 
rate coefficients from the minimal model of GP98. 



Radiative association of He and H + 

The rate given bvlGP98l was a fit to the rate coefficients of Roberge and 
Dalgarno Ird T 119821) "and Kimura et al. (KLDD9J). However, more accu- 
rate calculations are available from Jufek, Spirko and Kraemer (JSK95|) and 
Zygelman, Stancil and Dalgarno (ZSD98), which have been fitted by Stancil, 
Lepp and Dalgarno |SLD98^ for 10 K < T < 10 4 K, yielding 8.0 X 
10~ 20 ( r/300 Kr° 24 exp(-r/4000 K) cm 3 s - \ significantly below the rate of 
lGP98l (seeFig. B.l.a). 
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Stimulated radiative association of He and H + 

The rate of stimulated radiative association depends on the background radia- 
tion, which is assumed to be a blackbody like the CMB. It thus depends both 
on the g as tempe ratu re T and the radiation temperature T r . It has been com- 
puted by |JSK95l and lZSD98l and can be fitted for 10 K < T < 10 4 K with 
the formula 3.2 X 10 _20 r 1 - 8 (l. + 0.1(7/1 K) 204 )- 1 exp(-7 , /4000 K)(l. + 2 X 
10~ 4 (7V/1 K) 1 ') cm 3 s -1 . The rate is plotted for various radiation temperatures 
in Fig. B.l.b. 



Photodissociation of HeH + into He andYt 

This rate can be determined from detailed balance of the reverse re- 
action (radiative association). Using the rate of [JSK95 for radiati ve as- 
sociation yields a photodissociation rate slightly different from GP98: 
220(7V/1 K)°- 9 exp(-22740 K/T r ) s" 1 (c.f. Fig. B.2.a). 



Proton transfer: HeH + + H 



He + H+ 



IGP98I adopted a constant rate for this reaction, which w as based on a rate 
determination by Karpas, Anicich and Huntress (KAH79) using the ion cy- 
clotron resonance technique. A numerical integration of the new cross sec- 
tions of ILJB95 allows to give an improved rate with a slightly lower 
value, which can be fitted for 10 2 K < T < 10 5 K, yielding 0.69 X 
10" 9 (r/300 K)° 13 exp(-T/33100 K) cm 3 s -1 (see Fig. B.2.b). 7 
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Fig. B.l. New rate coefficients for helium chemistry, a) H + + 
He — > HeH + + y (rad. ass.). Dashed lines: Radiative association 
and inverse predissociation of lGP98l dotted line: ISLD98l b) H + + 
He + y — > HeH + + y. Stimulated association rate for black body 
radiation backgrounds with T r = 0, 500, 1000, 2000 and 5000 K. 



7 Note the discrepancy with lSLD98t who find a rate for the proton 
transfer reaction of HeH + which is larger by a factor of 1.5. 
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Fig. B.2. New rate coeffic ients fo r helium ch emistry, a) HeH + + 
y -» He+H + . dashed line: |CT9l dotted line: [SLD98l HeH + + 
H -hj He + H+. Soli d line: new fit, d otted lin e7 TsLD98l dashed 
line: lGP98l dot: lRDl(ll982h . triangles: OH1 



Table B.l. Collisional and radiative rates. T denotes the gas temperature in Kelvin, T eV the gas temperature in eV, T r the temperature of radiation in K. dex(x) = 10*. Acronyms: AAZN97: Abel, 
Anninos, Zhang, Norman (1997), CCDL07: Capitelli, Coppola, Diomede, Longo (2007), GJ07: Glover and Jappsen (2007), GP98: Galli and Palla (1998), GP02: Galli and Palla (2002), JSK95: 
Jufek, Spirko, Kraemer (1995), LJB95: Linder, Janev, Botero (1995), LSD02: Lepp, Stancil, Dalgarno (2002), MSM96: Martin, Schwarz, Mandy (1996), PSS83: Pequignot, Petitjean, Boisson C. 
(1991), SKHS04b: Savin, Krstic, Haiman, Stancil (2004, Erratum), SLD98: Stancil, Lepp, Dalgarno (1998), TT02: Trevisan and Tennyson (2002), ZDKL89: Zygelman, Dalgarno, Kimura, Lane 
(1989), ZSD98: Zygelman, Stancil, Dalgarno (1998) 



Reaction number Reaction 



Rate [cgs] 



Reference 
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2 

3 

4 

5 

6 

7 

8 

9 
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11 

12 
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14 
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16 
17 
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19 
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21 
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34 
35 
36 
37 
38 
39 
40 
41 
42 
43 



H+e-^ 
ir+H 



H+H+ 



H~+y 
-> H 2 +e 

->H+ +y 
Hj+H -> H 2 +H + 
2H+H -J H 2 +H 
H 2 +H+ -> Hj+H 
H 2 +e -> 2H+e 
H 2 +H -> 3H 
H~+e -> H +2e 
H~+H -> 2H+e 
H~+H + -> 2H 



exp(-T/16200) 



H~+H + 

Hj+e- 
H++H~ 
H 2 +e -> 
H~+y- 

Hj+r- 

H 2 +y — 

H 2 +y — 
D-+y - 
HD + +y 
HD + +y 
HD + +y 
HD+y- 
D + +e- 
D+H+ - 
D++H- 
D+H^ 
D+H 2 - 



Hj+e 
2H 

^H+H 2 
H+H~ 

• H+e 
H+H+ 
Hj+e 

• 2H + +e 
(H 2 )* -> 2H 

• D+e 
->D+H+ 

H+D+ 
-> H++D++e 
> HD + +e 
D +y 
>D++H 
*D+H+ 
HD+y 
. HD+H 



HD++H HD+H+ 
D + +H 2 -» HD+H+ 
HD+H -> D+H 2 

HD+H + — > D + +H 2 
D+H + -> HD++y 
D + +H^HD++y 
HD + +e -> D+H 
D+e -> D" + y 
D++D" -> 2D 
H++D" -> D+H 
H~+D -> H+D~ 
D~+H -> D+H~ 
D~+H -> HD+e 



1.4 X 10" 

1.5 X 10-' for T < 300 
4.0 X 10- 9 T- 017 for T > 300 

dex(-19.38- 1.523 log 10 (r) + 

6.4 x lO" 10 

5x lQ- 29 T- [ forT < 300 
see reference 

1.91 x lO-'r 136 exp(-53407.1/r) 
Fit to data of MSM96 
see reference 
see reference 

1.40 x 10- 7 (77300r 0487 exp(7729300) 



1.1181og 2 CO- 0.1269 log' CO) 



-9 T -0.35 



for T < 8000 



09 for T > 8000 
0.5 



6.9 x 10 _S T 

9.6 x io~ 7 r 

2.0 x io~ 7 r 

5x Ur 6 T-°- 5 for T> 100 
3.67 xlO 1 ^ 2 - 28 ^-^) 

1.1 x 10- l T? 13 exp(-8823/r r ) 
1.63 x 10 7 exp(-32400/7/,) 

2.9 x 10 2 7/,. L56 exp(-178500/7/ r ) 
9.0 x 10'T r L48 exp(-335000/r,) 
1.13 x 10 s r r 0369 exp(-140000/T r ) 
estimated by rate 16 
estimated by half of rate 17 
estimated by half of rate 17 
estimated by rate 19 
estimated by rate 1 8 
3.6 X 10-' 2 (77300r - 75 



37.1/T)-3.31 x 10- 17 r 148 
exp(-33.0/7/) + 2.03 x lo-'r -0 - 332 



2 X 10- 10 r 04u/ exp( 
2.06xl0- 10 r - 3% 
see reference 

1.69 x 10- 10 exp(-4680/r + 198800/7/ 2 ) for T > 200 
9.0 x 10- u exp(-3876/7") for T < 200 
6.4 x 10" 10 

1.0 x 10" 9 (0.417 + 0.846 log 10 (T) - 0.1371og 2 (r)) 
5.25 x 10- n exp(-4430/r + 173900/7/ 2 ) for T > 200 
3.2 x 10- u exp(-3624/r) for T < 200 

1.1 X 10- 9 exp(-488/7/) 

dex(-19.38- 1.523 log 10 (r) + 1.118 log 2 (r) - 0.1269 log 3 <T)) 



dex(-19.38 - 1.523 log 10 (r) + 1.11 

7.2 x io- 8 r-° 5 

3.0 x 10-' 6 (77300) a95 exp(-779320) 
1.96 x 10- 7 (r/300)" 0487 exp(T/29300) 
1.61 x 10- 7 (77300r 0487 exp(7729300) 
6.4 x 10- 9 (r/300) 041 

6.4 x 10- 9 (77300) 041 

1.5 x 10" 9 (r/300r° 1 



ilogt (D- 



0.1269 log)o(D) 



SLD98 

Savin (2002) 

Savin (2002) 

Dickinson (2005, 2008) 

GP02 

GP98 

SLD98 

GP02 

GP02 

GP98 

GP02 

GP98 

GP98 

SLD98 

SLD98 

LSD02 

LSD02 

SLD98 

SLD98 

SLD98 



Table B.l. Collisional and radiative rates in cgs units. T denotes the gas temperature in Kelvin, T eV the gas temperature in eV, T r the temperature of radiation in K. dex(x) = 10 v . Acronyms: AAZN97: 
Abel, Anninos, Zhang, Norman (1997), CCDL07: Capitelli, Coppola, Diomede, Longo (2007), GJ07: Glover and Jappsen (2007), GP98: Galli and Palla (1998), GP02: Galli and Palla (2002), JSK95: 
Jufek, Spirko, Kraemer (1995), LJB95: Linder, Janev, Botero (1995), LSD02: Lepp, Stancil, Dalgarno (2002), MSM96: Martin, Schwarz, Mandy (1996), PSS83: Pequignot, Petitjean, Boisson C. 
(1991), SKHS04b: Savin, Krstic, Haiman, Stancil (2004, Erratum), SLD98: Stancil, Lepp, Dalgarno (1998), TT02: Trevisan and Tennyson (2002), ZDKL89: Zygelman, Dalgarno, Kimura, Lane 
(1989), ZSD98: Zygelman, Stancil, Dalgarno (1998) 



Reaction number 
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Reference 


44 


D+I-T -» HD+e 


estimated by rate 42 


This work 


45 


H"+D + -> D+H 


1.61 x l(T 7 (77300r 0487 exp(7729300) 


LSD02 


46 


T-Tp+^-i-p i T-To+ -1-t/ 

JT.C TC * 1LC i y 


3 36 x W)~ 10 T~°- 5 tT /loom -02 ( i + it /m 6 i°' 7 )~ 


Ten M9991 


47 


nc \y nc t c 


see reference 
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T-Tp+4.p — itIp_i_T/ 

J. l\s it '1 It T Y 

He+y — » He + +e 


cpp rpfprpnrp 


AA7N97 


49 


see reference 


AAZN97 


50 


He+H + -> He + +H 


4.0 X 10" 37 r 4 " 74 for T > 10000 

1.26 X 10- 9 r-° 75 exp(-127500/r) for T < 10000 


GP98 
GJ07 


51 


He + +H -> He+H + 


1.25 X 10" 15 X (r/300) 25 


ZDKL89 


52 


He+H + — > HeH + +y, radiative association 


8.0 X 10- 20 (7;/300)" - 24 exp(-7;/4000) 


SLD98 


53 


He+H + +y — > HeH + +y, stimulated radiative association 


3.2 x 10- 20 T 1S /(1 + 0.1T 204 ) * exp(-r/4000)(l + 2 x 10- 4 7/,!') 


JSK95, ZSD98 


54 


He+Hj -> HeH + +H 


3.0 X 10-'°exp(-6717/7/) 


GP98 


55 


He + +H -> HeH + +y 


4.16 x 10- 16 r-° 37 exp(-7'/87600) 


SLD98 


56 


HeH + +H -> He+Hj 


0.69 x 10- 9 (77300) al3 exp(- 7733100) 


LJB95 


57 


HeH + +e -> He+H 


3.0 x io- 8 (r/300)-° 47 


SLD98 


58 


HeH + +y -» He+H + 


2207/ r a9 exp(-22740/7V) 


JSK95 


59 


HeH + +y -> He + +H 


7.8 X 10 3 T r L2 exp(-240000/r r ) 


GP98 



